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EFFECT OF COEFFICIENT CHANGES ON STABILITY OF LINEAR 


RETARDED SYSTEMS WITH CONSTANT TIME DELAYS 

L. Keith Barker 
Langley Research Center 


SUMMARY 

In previous studies, the effects of time delays on the stability of 
differential-difference equations of the retarded type (retarded systems) have 
been examined. After the stability condition of the retarded system is deter- 
mined for a fixed set of time delays, it may be of interest to determine the 
further change in the stability condition as a system coefficient is varied, 
with the time delays held fixed. This latter problem is examined in this paper 
by casting it in a form for which the x-decomposition method is applicable. The 
approach involves a rearrangement of the characteristic equation so that the 
coefficient plays a role similar to a time delay in the x-decomposition method. 
Then, the x-decomposition method is applied to obtain the changes in the stabil- 
ity condition as the coefficient is varied. 

The method is applied to a second-order differential equation with constant 
time delays in the velocity and displacement terms. Intervals of the coeffi- 
cient over which the system is stable and unstable are computed. 


INTRODUCTION 

A differential-difference equation of the retarded type (retarded system) 
is a differential equation in which the dependent variable and all but the high- 
est order derivative may contain time delays. (There is at least one delay.) 
Mathematical treatments of retarded systems can be found in references 1, 2 , 
and 3, for example. 

Some physical problems are modeled as a homogeneous retarded system, and 
the stability of the system is desired. Also, the stability of the homogeneous 
retarded system is related to the asymptotic stability of some nonlinear control 
problems involving time delays (ref. 4) and to the valid application of general- 
ized harmonic analyses to retarded systems (ref. 5). 

A homogeneous retarded system, which is linear and time- invariant , is asymp- 
totically stable if and only if all the roots of the characteristic equation 
have negative real parts (ref. 1). The characteristic equation is transcenden- 
tal and has an infinite number of roots. Therefore, it is not possible to com- 
pute all these roots to determine whether they have negative real parts. How- 
ever, it is possible to compute the number of roots with zero or positive real 
parts . 
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With zero delays, the characteristic equation reduces to a polynomial, 
which has only a finite number of roots. With sufficiently small positive time 
delays, the characteristic equation has essentially these same finite roots, 
plus an infinite number of roots with arbitrarily large negative real parts 
(refs. 6 and 7). Hence, with sufficiently small positive time delays, the sta- 
bility condition (stable or unstable) of the retarded system is the same as the 
stability condition of the system with zero delays, and the roots with nonnega- 
tive real parts are known. As the delays are increased from zero, the roots of 
the characteristic equation begin to move, generating root-locus curves. In 
order for the stability condition to change, a root-locus curve must intersect 
the imaginary axis. The x-decomposition method of reference 8 has been shown to 
be a convenient means for determining the directions in which the root-locus 
curves cross the imaginary axis. Hence, the change in the stability condition 
can be computed as the delays are varied, one at a time, from zero to their 
final desired values (ref. 7). 

Another stability problem occurs whenever a coefficient in the character^ 
istic equation of a retarded system must be varied. The purpose of this study 
is to show how the r-decomposition method may be used to solve this stability 
problem. The proposed technique involves a rearrangement of the characteristic 
equation so that the coefficient to be varied appears as a time delay. Then, 
the i -decomposition method is applied to obtain changes in the stability condi- 
tion as the coefficient is varied over a range of values. 


SYMBOLS 

A, A. N x N matrices of real constants 

f(s;x J1 ) function of s in equation (3) 

I N x N identity matrix 

Im( ) imaginary part of ( ) 

i imaginary unit, T 

l integer index (subscript) 

K system coefficient 

K m final desired value of K 

K* value of K at a point of intersection of a root-locus curve with 

the imaginary axis 

L( s) characteristic quasi-polynomial 

N dimension of x(t) in equation (1) 

N(K) number of roots of W(s) with positive real parts 
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X. + 


N( t ) number of roots of W 2 (s) with positive real parts 

Re( ) real part of ( ) 
s complex variable, a + ioo 

T number of time delays in system 

t time 

W(s) function of s in equation (7) 

W 2 (s) function of s in equation (6) 

x scalar function of time 

N x 1 state vector 
ot-j , ot 2 real numbers 

e small positive number 

o real part of complex number s 

t , t £ constant time delays 

t* value of x at a point of intersection of a root-locus curve with 

the imaginary axis 

a) imaginary part of complex number s 

u) m an upper bound on oj in L(s) = 0, where s = iw 

w* value of a) at a point of intersection of a root-locus curve with 

the imaginary axis 


ANALYSIS 
Retarded System 

A class of homogeneous dynamical systems, called retarded systems, can be 
described by the following equation: 

4* T 

x(t) = A x(t) + !> A^ x ( t— ) (1) 

£=1 

->■ 

where x(t) is an N x 1 vector, A and A^ are N x N constant matrices, 
and > 0 are constant time delays. 
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The characteristic equation associated with equation ( 1 ) is 


L(s) = det si - A - 


_T_ 

>_ 

A=1 


k SL e 


- t a s ' 


( 2 ) 


where s = o + iw is a complex variable. It is known (ref. 1) that x(t) -*• 0 
as t -*■ ® (asymptotically stable) if and only if all the roots of equation (2) 
have negative real parts (a < 0). 


The roots of equation (2) occur in complex-conjugate pairs, so that only 
roots in the upper half of the complex s-plane are considered herein; that is, 
u > 0. 


Change in Stability With Coefficient Changes 

It is assumed that the retarded system is known to be either stable or 
unstable for a fixed set of time delays □ 1, 2, . . ., T), and that the 

number of roots with positive or zero real parts has been determined. This 
information can be obtained by using the x-decomposition method of stability 
analysis (ref. 7). Now, suppose it is desired to know how the stability further 
changes as a system coefficient K is changed. Toward this end, the character- 
istic equation (2) is written as 

f(s;x A ) = Ks (3) 

where K is the system coefficient of interest, and f(s; Tjl ) contains the 
remaining part of the characteristic equation. The characteristic equation 
can be solved for K and the resulting equation multiplied by s to obtain 
equation (3). 

In order for the number of roots with positive real parts of equation (3) 
to change as K is continuously changed, a root-locus curve must intersect the 

imaginary axis. At an intersection point (s □ iw) , the following equations must 

be satisfied: 

Re[f(i(o;T £ )] = 0 (4) 

and 

K = - Imtfdujx.)] (5) 

to 16 

Equations (4) and (5) are the real and imaginary parts of equation (3) with 
s = iw. Thus, the points where the root-loeus curves of equation (3) cross the 
imaginary axis (eq. (4)) and the values of K (eq. (5)) for which these inter- 
section points occur can be calculated. The behavior of the root-locus curves 
at the intersection points determines the change in the number of roots with 

positive real parts. Specifically, it is desired to know in which directions the 

root-locus curves cross the imaginary axis or if a root-locus curve is tangent 



to the imaginary axis. These are determined later by using the x-decomposition 
method . 

The essence of the T-decomposition method of stability analysis is now pre- 
sented. Consider the function 

W 2 (s) = e TS (6) 

where W 2 (s) is a ratio of two polynomial equations in s with constant coeffi- 
cients, and x is a time delay. The following theorem is attributed to Lee and 
Hsu (ref. 8), where N(x) denotes the number of roots of equation (6) with posi- 
tive real parts. The theorem has been reworded for the present application. 

Theorem 1. Let s = ito* , where w* >0, be a purely imaginary root of equa- 
tion (6) with corresponding time delay x* ~ 0. Let 4> be the corresponding 
point of intersection of the testing function W 2 (iw) with the unit circle in 
equation (6). Then, 

(1) N ( x ) increases by 1 if the testing function enters the unit circle at 
$ (which occurs when x = x* and w = w*); 

(2) N(x) decreases by 1 as x increases across x* if the testing func- 
tion leaves the unit circle as oi increases across w*; and 

(3) N ( x) remains the same as x increases across x* if the testing 
function remains on the same side of the unit circle as w increases across 

0 )*. 


Lee and Hsu used the results of theorem 1 to examine the stability of a 
retarded dynamical system with one constant time delay x. The development of 
the theorem is of such generality, however, that the only restriction required 
on the testing function W 2 (s) in equation (6) is that it be analytic (or reg- 
ular) at the intersection points being considered. A significant feature of the 
theorem is that it is not necessary to be concerned with finding the first non- 
zero derivative of the real part of the root-locus curves with respect to K, or 
whether the intersection points are simple. 

Now, equation (3) is expressed in a form for which the x-decomposition 
method is applicable. Taking the exponential of both sides of equation (3) 
results in 


W(s) = e Ks 


(7) 


where 


f(s ; t ^ ) 

W(s) = e (8) 

Equation (7) is in the standard form for application of the x-decomposition 
method, where K is playing the role of a time delay; and a modified testing 
function has been introduced in equation (8) . 
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The following theorem results from the general proof of theorem 1 with 
equation (6) replaced by equation (7). 

Theorem 2. Let s = iio* be a purely imaginary root of equation (7) when 
K = K*. Moreover, let a-j < w* and ot 2 > w* be real numbers for which W(ia^) 
and W(ia 2 ) are defined and let w* be the only intersection point in the 
interval [a^c^]. Then for e an arbitrarily small positive number, 


(1) 

N(K*+e) = N(K») + 1 

if 

IW(ia 1 )| 

> 1 and | W( ic *2 ) I < 

1 

(2) 

N (K*+e) = N(K*) - 1 

if 

|W(ia 1 ) | 

< 1 and | W(ia 2 ) 1 > 

1 

(3) 

N(K*+e) = N (K*) if 

both 

1 W( ia 1 ) I 

1 and |W(ia 2 )l are 

greater than 1 , 


or if both are less than 1 . 

The values of K* at the intersection points are ordered by increasing magni- 
tudes to obtain the change in the stability as K increases to its final 
desired value, say K m . 

Geometrically, items (1), (2), and (3) of theorem 2 correspond, respec- 
tively, to the testing function W(iw) entering, leaving, and being tangent to 
the unit circle. The testing function is evaluated along the imaginary axis to 
infer the directions in which the root-locus curves cross the imaginary axis. 

Theorem 2 is applied as follows: Let s = iw* be a point of intersection 

of a root-locus curve of equation (3) with the imaginary axis when K = K*. 

Thus, to = co* and K = K* satisfy equations (4) and (5). This root-locus 
curve is common to both equations (3) and (7). In addition, there are no extra 
root-locus curves of equation (7), which are not common to equation (3), that 
have this intersection point when K = K*. If co = <o* is the only value of w 

in the interval [a^c^] for which the testing function W(ico) intersects the 

imaginary axis, then the change in the stability at this intersection point is 
determined by computing I W( i ) t and ) W ( i ot^ ) I . For example, from item (1) of 

theorem 2, if |W(ia^)| > 1 and |W(ia 2 )| < 1 , as indicated in figure 1, then 

the system gains exactly one root with positive real part as K increases 
across K*; that is, N(K*+e) = N(K*) + 1. 



Figure 1.- Hypothetical variation of the absolute value of the 
testing function, evaluated along the imaginary axis. 


The magnitude of w(iw) in figure 1 intersecting the horizontal dashed 
line at 1 corresponds to W(icj) intersecting the unit circle . Hence, in fig- 
ure 1, W(ioj) enters the unit circle as w increases across w* . 


Procedure 

To apply the stability procedure of this paper, w is incremented in the 
following equations: 

s □ iw (9) 


Im(K) 


Im 


s 


( 10 ) 


Re(K) 


Re 


f ( s;xp ; 
. 3 J 


(ID 


I f(s ; x £ } I 

I W( s) I = I e I (12) 

The coefficient K in equation (3) is assumed to be a real number. If 
Im(K) =0 in equation (10) for some value of w during the incrementation, say 
w* , then equation (3) is satisfied by s □ iw* when K in equation (11) has 
the real value denoted by K* . Theorem 2 is applied by examining tabulated val- 
ues of equation (12) as w increases across w* . 

Since w > 0 is incremented in equations (9) to (12), there is the ques- 
tion of when to terminate the iteration process. It is alway s possible to com- 
pute an upper bound on w for K ^ K m . The procedure is illustrated in a subse- 
quent application. It is important to choose the increment on w small enough 
to preclude overlooking any values of s = iw which satisfy the characteristic 
equation when K 1 K ffl . 


APPLICATION 

Consider the following equation from reference 6 (in the present notation): 


d 2 x dx, 

— (t) + 2K — (t-x ) 

dt 2 dt 1 


+ X(t-T 2 ) □ 0 


(13) 


where t ^ and t 2 are fixed time delays, and the real number K > 0 is a sys- 
tem coefficient. 
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The characteristic equation associated with equation (13) is 


— t 1 s -t 2 s 

s 2 + 2Kse + e = 0 04) 

With K = 0.5, t - j =0.8, and t 2 = 0.2, the solution of equation (13) is asymp- 
totically stable; that is, all roots of equation (14) have negative real parts 
(refs. 6 and 7) . 


Now, as K is varied, root-locus curves may possibly intersect the imagi- 
nary axis and change the stability of equation (13). To examine this possibil- 
ity, write equation (14) as 


Ks 



-t 2 s\ t i8 
+ e )e 


(15) 


which is in the form of equation (3) with 

. if 2 " T 2 s ) T 1 S 

f(s;T a ) = - -\s 2 + e /e 


( 16 ) 


Then, with equation (16), equations (9) to (12) along with the x-decomposition 
theorem are subsequently used. 


Iteration Bounds on w 

Maximum bound .- Setting s = itu in equation (14) gives 

-iWT. -itUTp 

-to 2 + 2Kicoe + e = 0 (17) 

Solving equation (17) for to 2 and taking the magnitude of both sides of the 
resulting equation gives 

I tol 2 1 2 |ku| + 1 (18) 

But, since u> > 0 and K > 0, 

to 2 £ 2K to + 1 < 2K m u +1 (19) 

where K < K ffl , and K m is a bound on K. 

When a) = 0, the right side of equation (19) exceeds the left side; how- 
ever, as w increases, the left side eventually dominates the right side. 

Thus, there exists a maximum value of to, say u> m , beyond which equation (19) 
is no longer satisfied. It is not difficult to see that this maximum value of 
to is the largest positive real root of the equation 

to 2 - 2K m to -1=0 (20) 
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which is 


m 


m 




m 


+ 1 


(21 ) 


Lower bound ,- One lower bound on a> is zero; however, a larger lower bound 
may be desired. In this case, equation (17) can be used to obtain the equation, 
or inequality, of magnitudes 


-ilDTp -iOJT- 

|u£| t e - 2Kiwe 


( 22 ) 


or. 


w 2 ^ 1 - 2Kio (23) 

Clearly, for w = 0, equation (23) is not satisfied. As to increases, there is 
a least value of u> for which equation (23) is satisfied.- This least value 
will be smaller for larger constant values of K. Thus, for K < K , the least 
value of oi which satisfies 

w 2 = 1 - 2K m <o (24) 

is taken as a lower bound on w. This lower bound, called , is the least 
positive real root of the equation 

to 2 + 2K m <0 -1=0 (25) 


which is 


“L - -K, * K* * ' 126) 

Hence, as K is increased to K m , all root-locus curves must intersect the 
imaginary axis in the interval 

W L - w = % (27) 


Specific Computations 

As mentioned previously, equation (13) is asymptotically stable when 
= 0.8, Tp = 0.2, and K = 0.5. It is desired to determine the change in 
stability as K is varied from 0.5 to 1 with the values of and Tp held 

fixed. In this case, eauation (16) becomes 

f ( s; t^) = - ^(s 2 + e -0.2s) e 0.8s (28) 

and equation (27) becomes 

0.41 < a) <2.5 (29) 
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With f(s;x a ) given by equation (28), equations (9) to (12) were iterated 
simultaneously on a digital computer by incrementing to in increments of 0.01 
over the interval of values in equation (29). 

Table I shows the regions of incremented values where Im(K) changes sign. 
The values of to where Im(K) = 0 and the corresponding values of K are 
given approximately in the first two columns of table II as to = to* and 
K = K*. The results in the third column of table II were determined by examin- 
ing |W(ito)| in table I in the neighborhood of to* and applying the 
x-decomposition theorem. As mentioned previously, the retarded system is known 
to be asymptotically stable for K = 0.5. Thus, for K between approximately 
K* = 0.16 and 0.63, the system will remain asymptotically stable. The system is 
unstable in the approximate intervals 0 5 K <0.16 and 0.63 < K = 1. To exam- 
ine values of K greater than K m = 1 , larger values of w m in equation (21) 
may be used. 


CONCLUDING REMARKS 

Once the stability condition (stable or unstable) of a retarded system is 
determined for a fixed set of time delays, it may be of interest to determine 
the further change in the stability condition as a system coefficient is varied, 
with the time delays held fixed. This latter problem is examined in this paper 
by casting it in a form for which the x-decomposition method is applicable. The 
approach involves a rearrangement of the characteristic equation so that the 
coefficient to be varied plays the role of a time delay. The x-decomposition 
method is applied then to obtain changes in the stability condition as the coef- 
ficient is varied over a range of values. 

The method has been applied to a second-order differential equation with 
time delays in the velocity and displacement terms. Ranges of values of a coef- 
ficient for which the system is stable and unstable were computed. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
January 13, 1977 
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TABLE I.- INCREMENTAL VALUES 


llllllllllllllllllllllll II 



Im(K) 

0.9887 


.0135 

.0087 

.0040 

-.0005 

-.0048 

-.0090 


-.0103 

-.0066 

-.0028 

.0012 

.0052 

.0093 


.5343 


|W(iw) I 
0.6667 


.9854 
.9905 
.9956 
1 .0005 
1 .0055 
1.0103 


1.0180 
1 .0115 
1 .0049 
-9979 
.9908 
-9835 


.2629 


TABLE II.- VALUES OF to*, K* , AND N(K*+e) 

0) = CO* 

1 .12 

1 .77 

12 NASA-Langley, 1977 L— 11298 
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